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Abstract 

We study the quark-antiquark potential of quenched SU(3) lat¬ 
tice gauge theory with the Liischer-Weisz action. After blocking the 
gauge fields with the recently proposed hypercubic transformation we 
compute the Sommer parameter, extract the lattice spacing a and set 
the scale at 6 different values of the gauge coupling in a range from 
a = 0.084 fm to 0.136 fm. 
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Recently a noticeable revival of the interest in improved gauge actions 
took place. It was observed that improved gauge actions make the numerical 
problems less severe when implementing chiral fermions. The underlying 
mechanism for the numerical improvement is a suppression of ultraviolet 
fluctuations of the gauge field. An example is the use of the Iwasaki and 
other improved gauge actions for domain wall fermions []'l] and also using 
the perfect gauge action || is an integral part of constructing the fixed 
point Dirac operator Q. Recently a systematic expansion of a solution of 
the Ginsparg-Wilson equation Q, the so-called chirally improved fermion, 
was proposed and implemented [|]. Also there it was found that using the 
improved gauge action is numerically advantageous for the Dirac operator. 
In particular the Liischer-Weisz action Q with coefficients from tadpole 
improved perturbation theory SI was used. Subsequently the instanton 
content of the QCD vacuum was studied for the Liischer-Weisz action in p]. 

In this letter we report on our results for the static potential and the 
lattice scale for the Liischer-Weisz action in order to make the use of this 
action easily accessible to the community. We furthermore test the recently 
proposed method of hypercubic blocking [|Io| which was found ]K), 11] to 
improve the statistical accuracy in the determination of the static potential 
by an order of magnitude. We analyze quenched ensembles at 6 different 
values of the gauge coupling and compute the Sommer parameter [12, [Tp 
and the lattice spacing a. This results in a precise determination of the 
lattice scale in a range between a = 0.084 fm to a = 0.136 fm and using a 
fit to our data even beyond this interval. Together with our results for the 
secondary couplings /?2 and /?3 this letter provides all ingredients necessary 
for using the Liischer-Weisz action at typical lattice spacings of state of the 
art simulations. 

In addition to the plaquette term of the Wilson gauge action, the Liischer- 
Weisz action includes a sum over all 2x1 rectangles and a sum over all 
parallelograms, i.e. all possible closed loops of length 6 along the edges of 
all 3-cubes. Explicitly the action reads 
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where (3\ is the principal parameter while /?2 and /3^ can be computed 
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from f3\ using one loop perturbation theory and tadpole improvement [Q], 
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The couplings P 2 , P 3 are determined self-consistently from uq and a for a 
given Pi. In Table |l] we list the values of the Pi used for our ensembles and 
our results for the expectation value of the plaquette no 4 = Re Tr(U p i)/3. 
The sample size at each value of Pi is 200 configurations on 16 4 lattices. 
The update was done with a mix of Metropolis and over-relaxation sweeps. 


Pi 

8.00 

8.10 

8.20 

8.30 

8.45 

8.60 

V 

P 2 

p3 

0.62107(3) 

-0.54574 

-0.05252 

0.62894(3) 

-0.54745 

-0.05120 

0.63599(3) 

-0.54998 

-0.05020 

0.64252(3) 

-0.55332 

-0.04953 

0.65176(3) 

-0.55773 

-0.04829 

0.66018(3) 

-0.56345 

-0.04755 


Table 1: Parameters for the Liischer-Weisz action. We list the values of the 
Pi and the expectation value of the plaquette uq 4 = Re Tr(f7 p /)/3. 


Before measuring the potential we applied the hypercubic blocking trans¬ 
formation proposed in (lCj. Hypercubic blocking mixes only gauge links 
from the hypercubes attached to the target link and has less impact on the 
short distance properties of the gauge fields than previously used smearing 
methods. The hypercubic blocking transformation proceeds in three steps 


Vi,fi-,vp ~ 'PsUi 3) 


(1 CK3)Ui tfl + 2 X/ Ui )V Ui+r),p,Uj +l ^ 
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In the first step intermediate fields Vi^- )V p are created from the thin-link 
variables Ui )P (indices i run over all sites of the lattice and fi, zz, p and p over 
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the four directions). In the second step the intermediate fields Vi tfl - U p 
blocked into a second set of intermediate fields which in the third step 
are transformed into the final fields 14 i/i . The restrictions on the indices fj,, v 
and p implemented in the sums in Eqs. (Q) ensure that contains only 
contributions from the hypercubes attached to the link (i,p). By Vgu( 3 ) we 
denote the projection of the sums back to elements of SU(3). The parameters 
aq, «2 and «3 determine the admixture of staples in each step of the blocking 
process. These parameters were optimized |l(J to minimize the fluctuations 
of the plaquette. Their values are given by oq = 0.75, «2 = 0.6 and 03 = 0.3. 

We measured the static potential on the smeared configurations using 
planar Wilson loops W (r, t) of size r x t with both t and r ranging from 1 
to 10. We fitted the expectation values of the Wilson loops to a sum of two 
exponentials c\ exp (-V(r)t) + C 2 exp(— E't) in a range of t = 2, 3 ... 9. The 
second exponential takes into account the contribution from excited states 
E' > V{r) and from the first term we directly obtain the potential V(r) 
for two static sources at distance r. As a cross check we also computed for 
some of the ensembles the potential for the raw, unblocked configurations. 
We find that the results are compatible within error bars but the statistical 
fluctuations, in particular at larger values of r and t, are much more severe 
for the raw configurations. 
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Figure 1: The static potential for the (5\ = 8.60 ensemble. The diamonds 
represent the raw data from the two parameter fit on the Wilson loops 
and the bursts are the values corrected for the short distance effect of the 
hypercubic blocking. The full curve is the parameterization (j|). 
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In Fig. |l] we show our results for the static potential for the (3\ = 8.60 
ensemble. The smooth curve is the standard infrared parameterization for 
the continuum potential, 

V(r) = C — A/r + or , (5) 

with constants C, A and a. The diamonds are the values for the potential 
obtained from the Wilson loops. The error bars are smaller than the sym¬ 
bols. For small distances, one finds a noticeable deviation from the Coulomb 
behavior —A/r. This deviation is an effect of the hypercubic blocking. How¬ 
ever, this effect can be computed perturbatively and the obtained deviation 
from the Coulomb potential is used to introduce a fourth fit parameter in 
the potential fit |ll| . Subtracting this perturbative part gives the corrected 
data which we represent by bursts. From the parameters A and a we com¬ 
puted the Sommer parameter ro jf^] in lattice units a and assuming ro = 0.5 
fm we extracted the lattice spacing a. The result is very stable under vari¬ 
ation of the fit range for the potential. Even the r/a = 1 measurement 
can be included. The final result is the weighted average over all fit ranges 
[ar min , ar max ] with r min G {1,2, 3} and r max G {7,8, 9}. We give the results 
for the lattice spacing and the Sommer parameter in Table 


Pi 

8.00 

8.10 

8.20 

8.30 

8.45 

8.60 

ro/a 
a [fm] 

3.688(37) 

0.136(1) 

4.015(34) 

0.125(1) 

4.362(41) 

0.115(1) 

4.741(49) 

0.105(1) 

5.289(66) 

0.095(1) 

5.967(70) 

0.084(1) 


Table 2: Results for the Sommer parameter ro in lattice units and the 
corresponding values for the lattice spacing a when ro is assumed as 0.5fm. 


Without the improved fit ro can still be determined using (]||). In this 
case, the lower limit of the fit must be chosen so that the region where the 
HYP smearing distorts the potential is excluded, as can bee seen from Figure 
Q, i.e. r m i n G {2,3}. The results of an analysis based on the three parameter 
fit (|5) are in perfect agreement with those given in the above table (for e.g. 
pi = 8.30 one obtains ro = 4.732(43)). Thus the only benefit from the 
improved fit is that no data points have to be excluded. This shows that 
the effect of the HYP smearing on the short-distance static quark potential 
is under good perturbative control. 
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In order to make the Sommer parameter and the lattice spacing available 
also for other values of [3\ we fit our data to a functional form based on the 
/3-function as proposed in [|R| . We find 

ln(r 0 /a) = 1.55354 + 0.79840 {Pi - 8.3) - 0.09533 {Pi - 8.3) 2 . (6) 

In Fig. we compare our numerical data for r$/a and a (again assuming 
?’o = 0.5 fm) to the curve (§). It is obvious that the data are well described by 
our parameterization. Furthermore when extending the plot range to values 
of Pi as small as Pi = 6.8 we find that our results are in good agreement 
with the data computed for very coarse lattices in i.e. a = 0.24 fm at 
Pi = 7.4, a = 0.33 fm at Pi = 7.1 and a = 0.40 fm at Pi = 6.8. 




Figure 2: Results for the Sommer parameter in lattice units (left-hand side 
plot) and the lattice spacing in fermi (right-hand side plot) as a function of 
Pi. The full curves are the interpolations of the data with the function (j6j). 
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Figure 3: Superposition of the potentials for all values of B\. The irrelevant 
constant C in Formula (j^) is set to zero. 
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Finally in Fig. ^ we show a common plot of our results for the static 
potential at all values of f3\ we analyzed. We set the irrelevant overall con¬ 
stant C to zero for all (3\. It is obvious that the data from different lattice 
spacings are in perfect agreement and the discretization errors are hardly 
noticeable for the Liischer-Weisz action. 
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